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Abstract 

This paper is devoted to study the Bianchi type III model in the 
presence of anisotropic fluid in f{R) gravity. Exponential and power- 
law volumetric expansions are used to obtain exact solutions of the 
field equations. We discuss the physical behavior of the solutions and 
anisotropy behavior of the fluid, the expansion parameter and the 
model in future evolution of the universe. 
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1 Introduction 

Hubble diagram of type la Supernovae (SNela) measured by the Super- 
novae Cosmology Project [Ij and High- 2; Team, up to the redshift z ~ 1 
[2], has been the first piece of evidence that current universe is undergoing a 
phase of accelerated expansion. The other Balloon-born experiments such as 
BOOMERanG [5] and MAXIMA [1] have detected the anisotropy spectrum 
of the Cosmic Microwave Background (CMB) radiations representing that 
universe is spatially flat. These data indicate that present universe is domi- 
nated by an un-clustered fluid with large negative pressure called dark energy 
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(DE) which causes expansion. The above picture has been further strength- 
ened by the current measurements of CMB spectrum obtained by Wilkinson 
Microwave Anisotropy Probe (WMAP) experiment |5],[6] and by extension 
of the Hubble SNela diagram to redshifts higher than one |7]. There are 
various attempts to construct the acceptable dark energy models in different 
directions. For example, traditional cosmological constant, quintessence or 
phantom models, dark fluid with complicated equation of state. String or 
M-theory, higher dimensions, brane-world models, etc. In addition to these 
attempts, the f{R) theory of gravity has also been helpful in describing the 
evolution of the universe. 

The recent enthusiasm in f{R) theory is caused by its success as the grav- 
itational alternative for DE |S]. Although, the cosmological constant. A, is 
the simplest explanation of the DE and a best fit to some of the astrophysical 
data |5j. However, the A-Cold Dark Matter (ACDM) fails in explaining why 
the value of A is so tiny (120 orders of smaller magnitude) as compared to 
the vacuum energy predicted by particle physics (the coincidence problem). 
In order to solve this problem, people replaced the cosmological constant 
with a scalar field referred to quintessence or phantom models [H]. In order 
to explain DE of the current universe with the help of the above models, one 
can use the effective equation of state parameter (EoS) u. When w = — 1, 
the universe passes through ACDM epoch, li w < —1 then we live in the 
phantom-dominated universe and for w > —1, the quintessence dark era oc- 
curs. It is mentioned here that energy conditions are violated in all these 
cases. Also, the phantom phase ends at finite-time singularity in the future 
[lOj while quintessence universe may end up at more general singularity. This 
is also true for the f{R) gravity which leads to DE with the corresponding 
effective u (quintessence or phantom) [S]. 

The f{R) theory of gravity quite naturally describes the transition from 
deceleration to acceleration in such an evolving universe [11]. This theory is 
very useful in high energy physics, for example, to solve hierarchy or gravity- 
GUTs unification problems [12]. A large body of papers is available in the 
literature [13]- [U] addressing the well-known issues of gravitational stabil- 
ity [ini, Newtonian limit [T7|, singularity problems [IB], solar system test 
[19], etc, in the context of f{R) gravity. The exact solution of the field 
equations in metric f{R) theory has been obtained in all the three symmet- 
ric versions. Spherically symmetric solutions have been investigated for both 
vacuum and non- vacuum cases by different authors [20]. Azadi et al. [21] and 
Momeni [22j have studied vacuum cylindrically symmetric solutions. Sharif 
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and Shamir [23] have explored plane symmetric solutions both for vacuum 
and non- vacuum cases. 

It is found that some large-angle anomalies appear in CMB radiations 
which violate the statistical isotropy of the universe [21]. Plane Bianchi 
models (which are homogeneous but not necessarily isotropic) seem to be the 
most promising explanation of these anomalies. Jaffe et al. [2S] investigated 
that removing a Bianchi component from the WMAP data can account for 
several large- angle anomalies leaving the universe to be isotropic. Thus the 
universe may have achieved a slight anisotropic geometry in cosmological 
models regardless of the inflation. Further, these models can be classified 
according to whether anisotopy occurs at an early stage or at later times of 
the universe. The models for the early stage can be modified in a way to 
end inflation with a slight anisotropic geometry [26]. For the latter class, 
the isotropy of the universe, achieved during inflation, can be distorted by 
modifying DE [27] . 

The universe acceleration provides information about the major part of 
the universe which has large negative pressure but without telling anything 
about the number of cosmic fluids in the universe. This may be explained 
by considering the accelerating expansion with a single fluid and an equation 
of state acting like DE. The main benefit of this approach is that a suitable 
equation of state can be obtained and observational data can be fitted. In 
General Relativity (GR), people have worked by choosing anisotropic fluid 
with anisotropic EoS. Akarsu and Kilinc [28] have studied the Bianchi type III 
model in the presence of single imperfect fluid with dynamical anisotropic 
EoS parameter. They have concluded that anisotropy of the DE do not 
always promote anisotropy of the expansion. Sharif and Zubair [29] have 
investigated the Bianchi type VIq cosmological models in the presence of 
electromagnetic field and anisotropic dark energy. They have examined the 
effects of electromagnetic field on the dynamics of the universe. In a recent 
work [30], we have obtained solutions of the Bianchi type VIq universe in 
f{R) gravity. 

Here, our objective is to find exact solutions of the Bianchi type III model 
to discuss future evolution of the universe in the metric f{R) gravity. We 
take anisotropic fluid to represent DE. The paper is organized as follows. 
In section 2, we present some features of the Bianchi type III model and 
define dynamical quantities describing the evolution of the universe. Section 
3 provides the field equations and anisotropic parameter of the expansion. 
We obtain possible solutions of the field equation and discuss the physical 
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behavior of these solutions in section 4. In the last section 5, we summarize 
the results obtained. 



2 The Model 

The homogeneous and anisotropic Bianchi type III spacetime is described by 
the line element 

ds^ = df - A\t)dx^ - e-^''''B\t)dy^ - C\t)dz'', (2.1) 

where the scale factors A, B and C are only functions of cosmic time t and 
a 7^ is a constant. 

The energy-momentum tensor for anisotropic fluid is given as 

= diag[p, -p^, -py, -p^] = diag[l, -u^, -Uy, -uj^]p, (2.2) 

where p is the density of the fluid while Px, Py and Pz are pressures and uux, 
oOy and Uz are directional EoS parameters on the x^y and z axes respectively. 
The DE schematically characterized by this parameter. The deviation from 
isotropy may be obtained by setting 

uix — i^-i ujy — ui -\- 5 and = a? + 7, 

where ui is the deviation-free EoS parameter and 5 and 7 are the deviations 
from u on y and z axes respectively. In this case, the energy-momentum 
tensor becomes 

r; = diag[l, -uj, -{u + 5), -{u + j)]p. (2.3) 

The skewness parameters S and 7 may be constant or functions of cosmic 
time. The average scale factor and volume of the universe for this model will 
be 

a = (ABC)^/^, V = a^ = ABC. (2.4) 

The mean Hubble parameter H, deceleration parameter q, expansion scalar 
G and the directional Hubble parameters in the x, y and z directions turn 
out to be 
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Here the deceleration parameter q measures the rate of expansion of the 
universe. The sign of q indicates the state of expanding universe. If g < 
or g > respectively, then it represents inflation or deflation of the universe 
while g = shows expansion with constant velocity. 

To examine whether expansion of the universe is anisotropic or not, we 
define anisotropic expansion parameter as 

^ = ^E(^^)^ (2-7) 

1=1 

where Hi {i = 1,2,3) are the directional Hubble parameters in the direc- 
tions of X, y and z axes respectively. If A = 0, then the universe expands 
isotropically. Further, any anisotropic model of the universe with diagonal 
energy-momentum tensor approaches to isotropy if A — > 0, K — > +oo and 
p > as t ^ +00 ([28], [31]). 



3 The Field Equations in f{R) Gravity 

The Lagrangian in f{R) gravity is given as 

. f{R) 



2k 



+ Lm, (3.1) 



where f{R) is a function of the Ricci scalar and Lm describes all kinds of 
matter including non-relativistic (cold) dark matter. Assuming variation 
with respect to the metric tensor, g^^, we have the following fourth order 
partial differential equations 

F{R)R^, - ^f{R)g^u - V^V.F(i?) + g^u^FiR) = kT^,, (3.2) 

where F{R) = df{R)/dR, □ = V^V^ with representing the covariant 
derivative and k(= = 1) is the coupling constant in gravitational units. 
The trace of Eq. (13. 2 p yields 

F{R)R-2f{R) + 3nF{R) =T. (3.3) 

This equation is helpful to solve the field equations and also to express f{R) 
in terms of its derivative, i.e., 

m = 'T^F(R)R^3DF(R)^ ^^^^ 
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The scalar curvature for Bianchi type III model is given by 



R = -2 



A B C AB AC BC 

A^B^C~A^^AB^AC^^ 



(3.5) 



The corresponding field equations become 



,A B C,^ ,A B C,- 



,A AB 
'A~A^^AB 



B 
B 



a 



AB 
A^^AB 



F 

.A 



B 
B 



c 



= -p, (3.6) 
-wp, (3.7) 



,C AC EC 1 ^ 



A B' 



a{- - —)F 



'A B 



-{u + 6)p, (3.8) 
-(a; + 7)p, (3.9) 
0. (3.10) 



The solution of Eq.f l3.10p yields 



B = ciA, 



(3.11) 



where Ci is a constant of integration. Using Eq.f l3.11|) and then subtracting 
Eq. (l3.8p from Eq.( ]3.9p . we obtain 5 = 0. This indicates that directional EoS 
parameters Ux, (jJy along x and y axes become equal and hence pressures. 
Consequently, the field equations turn out to be 



,A 
'A 
,C_ 

c 



A^ AC,^ 1 ,A C . 



(3.12) 
—up, (3.13) 
(3.14) 



Subtracting Eq. (13.131) from Eq. fl3.14p and integrating the resulting equation, 
we have 

A C c 1 f ,a^F , , , , 

+ — / i^ + -fp)Vdt, (3.15) 
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where c is a positive constant of integration. In terms of directional Hubble 
parameters, the above equation can be written as 

Using Eqs.daH]) and (EZD, it follows that 

A = ^[c + / {^ + lp)VdtfV-'F~\ (3.17) 

If we take 7 = and F{R) = 1, the anisotropy parameter of expansion 
reduces to GR for an isotropic fluid 

The integral in Eq. fl3.17p vanishes for the following value of 7 [28l [29], i.e., 

7 = -^. (3.19) 

Consequently, the energy-momentum tensor, anisotropy parameter and Eq. fl3.16p 
will take the form 

= diag[l,-u,-u,-u + ^^]p, (3.20) 

A = l^,V-^F'^ (3.21) 
H.-H, = (3.22) 

It is interesting to mention here that these results will reduce to GR [28J for 
FiR) = 1. 



4 Solutions of the Field Equations 

In this section we would obtain exact solutions of the Bianchi type III model 
in the presence of anisotropic fluid. Using Eq. fl3.19p in Eqs. fl3.12p - (I3.14I) 
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along with Eq.f l3.3p . we have 

(4.1) 

+ + IC^^ + 2^^^^ - F - (- + -)F = -(a; + 7)^, (4.2) 

C 2AC 1 9/i 

+ -j^)F + -f{R) -F- {—)F =-{uj + 7)p, (4.3) 

,2i C 2^(7,^ 

+3(— + ^)F=(l-3a;-7)p. (4.4) 

This is a system of four differential equations in five unknowns. In order to 
obtain solution, we make use of two different volumetric expansion laws (in 
next subsections) to get a system of five equations with five unknowns. 

4.1 Bianchi Type III Model for Exponential Expansion 

We consider the following volumetric exponential expansion law 

V = C2e^^\ (4.5) 

where C2 and k are positive constants. Using the value of V in Eq. (13.221) . we 
have 

H..-H, = — . (4.6) 

C2F 

Solving the system of equations (I4.ip -( l44|) and Eq.f l4.6p . we obtain the scale 
factors as follows 

A = 036*^*+^^'^''*, (4.7) 

-3kt 

B = cicze''^~2 J (4.8) 

C = c^e'''-^ ^ (4.9) 

where C3 and C4 are positive constants of integration. To solve the integral 
part in the above equations, we assume a relation between F and V as 
FaV""/^ [23], which gives 

F = C5V"'/^ F = c^e'^^\ (4.10) 
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where C5 is a positive proportionality constant and m is a non-zero arbitrary 
constant. Using Eq. fl4.10p in Eqs.f l4.7p - fl4.9p . the scale factors become 

A = c^e ^^a'^s ^(3+™) ^ (4.11) 

B = cicge sc^c^WT^" , (4.12) 

C = C4e ^+3c2C5M3+m)^ . (4.13) 

This is the first solution of Bianchi type III model with exponential volu- 
metric expansion. When m > —3, the scale factors admit constant values 
at early times of the universe, after that start increasing with the increase 
in cosmic time without showing any type of initial singularity and finally 
diverges to 00 for t ^ 00. This shows that at the initial epoch, the universe 
starts with zero volume and expands exponentially approaching to infinite 
volume. However, for m < —3, the scale factors A, B increase with time 
while C approaches to zero as t — )■ 00. Also, the model with m = —3 rep- 
resents the universe at early stage. Moreover, the expansion scalar for these 
scale factors exhibits the constant value, i.e., = 3A; which shows uniform 
exponential expansion from t = to t = 00. 

The mean, directional and deceleration parameters turn out to be 

= n=k + ^e-^'^"'^'\ = k - ^e-^'^"^^'*, 

OC2C5 <JC2C5 

H = k, q = -l. (4.14) 

This shows that the mean Hubble parameter is constant whereas others are 
dynamical. As time approaches from zero to infinity (for m > —3), the 
directional Hubble parameters start reducing towards the constant value of 
H and becomes equal as t — > 00. For m < —3, parameters along x and 
y axes will increase from the mean Hubble parameter by a constant factor 
g^j^ while parameter along z-axis decreases by twice the same factor. Also, 
the deceleration parameter appears with negative sign which shows acceler- 
ating expansion of the universe as we can expect for exponential volumetric 
expansion. 

The anisotropy parameter of the expansion takes the form 

A = -(-^)V2(-+3)^*. (4.15) 
9 kc2C^ 

Here one can observe that at t = 0, the anisotropy parameter measures a 
constant value while it vanishes for m > —3 at infinite time of the universe. 
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This indicates that the universe expands isotropically at later times without 
taking any effect from anisotropy of the fluid. However, for m < —3, the 
anisotropy in the expansion will increase with time. 

Inserting the scale factors, Eqs.f l4.lip -f l4.13p . in Eqs. fl4T]) -f lO|) and use 
Eq.f l4.10|) . we obtain the following energy density of anisotropic dark fluid 

p = ^^[-3m2A;2 + 2mA;2 + -I^e-(3+-)'=*- -^6-2(3+-)'=* 

3 3 C2C5 {C2C5y 

- — 3C2C5 fc(3+m) ]_ (4.16) 

We see that the matter density is constant at early stage of the universe 
[t = 0) and shows monotonic behavior in the evolving cosmic times. This 
behavior of matter density deviates from GR where it becomes constant 
in exponential expansion. Here, for positive values of m, it will positively 
increase with time and diverges to oo, showing future time singularity. For 
—3 < m < 0, it starts decreasing with cosmic time and eventually approaches 
to zero as t — )■ oo. It is mentioned here that isotropic conditions of the model 
may not be satisfied and hence leading to the conclusion that Bianchi type 
III model remains anisotropic in the later times whereas in GR it becomes 
isotropic. The anisotropy of the model is maintained not due to the presence 
of anisotropic fluid but due to modification in gravity. It is worth mentioning 
here that the results correspond to GR results for F{R) = 1. 

The anisotropic EoS parameter obtained from the field equations is given 

by 

2 — (3 + m)fct 

u = [3m(fcc2C3C5)^ - 2mfccc2C5e-(^+'")'=* + 9a^clcle~^'"''^^^^^^]/ 
[3m{kc2C3C5)\-3m + 2) + 2mkcc2C5e-^^+"'^^' - 60^0^6-2(3+'")'=* 

2c g — (3+m)fci 

-Qa^cjcle' ^-'^2^5 Ma+m) ]. (4.I7) 

When t — )■ 00 this gives a constant quantity, 2-3m ' "W'hich can further be 
characterized into phantom and quintessence regions for m > lor— 3<m< 
1 respectively. For m = 1, the EoS parameter is w = — 1, corresponds to the 
vacuum energy which is mathematically equivalent to cosmological constant 
while m < —3 does not yield fruitful results. We see that p increases rapidly 
in the phantom region where a; < — 1. It is mentioned here that phantom 
regime support the observational evidence of a recent supernova data |33j - 
[35] . However, it decreases in the quintessence region representing relatively 
slow expansion. 
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Using Eqs.f l4.1(]]) . fl4.1ip and f l4.16p in Eq.f l3.19p . the skewness parameter 
along 2;-axis becomes 

[3m(A;c2C3C5)'(-3m + 2) + 2mkcc2C5e-^'^+"'^''' - 6c^cle-^^'^+"'^''' 
-Ga^CgCge" ^3+™) ]. (4.18) 

We have aheady found 5 = and we see that 7 also approaches to zero 
as t — > 00. Thus the anisotropy of the fluid is completely removed in the 
future evolution of the universe. The scalar curvature and f{R) function for 
exponential model becomes 

9 C2C5 3 C2C5 C3 

f{R) = ^^^[(SQm^e + Mme + 9R)clclcl - 4mA;Ac2C5e-'(3+™)^* 

o\ — (3 + m)fct 

+6A2c^e-2(=^+'")'=* + 24«2^2^2g-'*+^^T^]. (4.19) 

This becomes constant in the quintessence region when t — )■ 00 while for 
m < —3, it diverges giving the same behavior as the energy density. Since 
the energy density decreases in the quintessence region, thus the universe 
starts from inflation driven by the anisotropic energy density at the early 
stage where curvature is very large. The scalar curvature also reduces with 
the passage of time. After that time, the matter density or radiations become 
small and the curvature becomes constant. 



4.2 Bianchi Type III Model for Power-law Expansion 

Here we solve the field equations by assuming power-law volumetric expan- 
sion given as 

V = C2t^", (4.20) 

where C2 and n are positive constants. For n > 1, Bianchi models exhibit ac- 
celerating volumetric expansion. When n = 1, the models represent volumet- 
ric expansion with constant velocity while for n < 1 these show decelerating 
volumetric expansion. Adopting the same procedure as in the exponential 
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model, we obtain the scale factors 



l — Sn—mn 



A = Care 3-2-, (1-3— ™«)' ^ (4.21) 



l—3n—mn 



B = ciC3re'^''2''o(i-3-'"")' , (4.22) 

—'2c ,X — 2n—mn 

C = C4re3<=2'=5(l-3n-mn)^ _ (423) 

We discuss the behavior of these scale factors in the future evolution by 
taking very large values of time. It is found that for n < the scale 

factors A and B increase with time while C approaches to zero provided 
that m is always be greater than —3 because n is assumed to be a positive 
number. This indicates that the universe is expanding along x and y axes 
with deceleration whereas there is no expansion along 2;- axis. Similarly, for 
n > 3:^^, the behavior of scale factors interchanged, however, the range 
— 3 < m < — 2 lies in the accelerating phase and — 2 < m comprises the 
decelerating epoch. 

The Hubble parameters, scalar expansion and anisotropy parameter of 
expansion in this case become 

H = ^, H^ = ^ + '''''"'\ H^^^-^t-^^^-^-, (4.24) 
t t 3C2C5 t 3C2C5 

e = — , A = -(^— )V(^-3"— (4.25) 
t 9 nc2C5^ 

Wc sec that the Hubble parameters and scalar expansion are extremely large 
at the origin of the universe and start decreasing monotonically with the 
passage of time and possibly will take zero value in the future (provided that 
m > —3). This shows that, at earlier times of the universe (just after the 
big bang), the expansion is much faster but slows down for later time of the 
universe. The anisotropy parameter approaches to zero for n > leaving 
isotropic expansion in the future while for n < its behavior is switched. 
With the same manipulations as in the previous subsection, the unknown 
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functions in the field equations turn out to be 

3 3 C2C5 

0„2 9rV^ —2c tl-Sn-mn 

_(i:^)2^-2(3+m)n _ f^rL^-angsS^i (i_3n-^n) 1 ^ ^4 26) 

C2C5 c| 



/[3(c2C3C5)2(-37n7l(m7l - 1) + 277171^ + 671)^^ + 2771710020503^" 

— 2c jl — 3n — rn.n 

-6o'o^t-'(3+"')" - 6a2o^o^t-2"e3^2C5(i-3n-™n)]^ (4.27) 

— 2c jl — 3n — rn.n 

7 = [-9a2o205t-2"e3^(i-3"-'"")]/[2(o20305)2(-3m7z(m7z - 1) + 277171^ 

+67l)r2 + 277l7l0020502r(3+-)"-l - 6020^^2(3+-)- 



1—Zn—mn 



-2c t 



1— 3n— mn 



-6Q;''o205t~''"e3=2C5 (i-3n-mn)]_ (4 28) 

In the accelerating models, the energy density will increase with times while 
for the decelerating models it will decrease. The dark energy EoS parameter 
for the accelerating models gives co — —1 + 6-3"^("^"-i)+3"tn latter times. 

^ 6— 3m(mn— l)+2mn 

This indicates that the universe will pass through quintessence epoch with 
appropriate values of tti > —3. The skewness deviation parameter along pr- 
axis vanish in the future evolution. Thus DE isotropous in the far future. 
The scalar curvature is 

R = (671- 1071^)^2 --4!^r2(3+-)" + i^r(^+-)" 
^ ^ 9 oioi 30205 

2q/2 ^l—3n—mn 
+ __i-2ng3c2C5 (l-3n-mn)^ (4.29) 

^3 

f{R) = [(477i7i(7n7z — 1) H mn^ + 671)^-^ + R 

_|_ ^-2ng3c2C5 (i-Sn-mn) _| ^-2(S+m)n 

3 o| 3o3 902O5 

(4.30) 

The scalar curvature approaches to zero for the passing times similar to p. 
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5 Summary and Outlook 



The aim of this paper is to discuss expansion of the universe due to anisotropic 
DE in f{R) gravity. For this purpose, we have found exact solutions of 
Bianchi type III model by assuming exponential and power-law volumetric 
expansion to cover all types of the expansion histories. These models repre- 
sent an accelerated expansion of the universe with V ^ oo as t ^ oc and 
supports the observations of the Type la supernova [H [2] and WMAP data 
[5[ |6]. The whole discussion is made in terms of F{R) which is assumed to 
have direct proportion with volume of the universe. The results of the paper 
can be summarized as follows: 

• The solution of the field equations for Bianchi type III model yields that 
the scale factors A and B are equal and the deviation from isotropy 
along y-axis is zero (i.e., 6 = 0). This shows that pressure of DE along 
X and y axes is same. 

• The physical behavior of the dynamical quantities depends on the val- 
ues of m. In case of exponential model, the dynamical parameters show 
different behavior for m ^ —3 while in power-law model we can only 
discuss m > —3 for later times to keep n positive. However, here we 
have further two types of models for n < or n > 

• The scalar expansion is constant for the first solution indicating that 
the universe expands homogeneously whereas in the second solution 
it shows that the expansion rate is faster at the beginning but slows 
down with the passage of time. The anisotropy parameter of expansion 
at t = measures a constant value while it vanishes for m > —3 at 
infinite time of the universe. This shows that the universe expands 
isotropically at later times. It is interesting to mention here that our 
solution can be used to understand cosmological phase transitions. It 
may help to understand the isotropization mechanism which could be 
responsible for the passage from a possible prior anisotropic phase to 
the present isotropic era we live in. 

• Behavior of EoS parameter and p exhibit that the universe will pass 
through quintessence region and phantom regions for m > 1 or —3 < 
m < 1 respectively. It is worth mentioning here that solution in the 
phantom regime {u < —1) supports the observational evidence of a 
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recent supernova data [23]- [SS]- Notice that the energy density in- 
creases/decreases monotonically with the passage of time depending 
upon the value of m in f{R) gravity. However, it becomes constant in 
GR at later times. The skewness parameter along z-axis will vanish 
thus acquiring isotropy in the fluid at infinite times. 

• The Bianchi type III model remains anisotropic even in the later times 
which is different from GR [28j where these become isotropic. 

• The scalar curvature becomes constant for m > —3 and diverges for 
m < —3 in exponential models while it becomes zero in the power-law 
models for later times. 
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